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ABSTRACT 

We present broadband (gamma-ray, X-ray, near-infrared, optical, and radio) observations of the 
gamma-ray burst (GRB) 090709A and its afterglow in an effort to ascertain the origin of this high- 
energy transient. Previous analyses suggested that GRB 090709A exhibited quasi-periodic oscillations 
with a period of 8.06 s, a trait unknown in long-duration GRBs but typical of flares from soft gamma- 
ray repeaters. When properly accounting for the underlying shape of the power-density spectrum of 
GRB090709A, we find no conclusive (> 3c) evidence for the reported periodicity. In conjunction 
with the location of the transient (far from the Galactic plane and absent any nearby host galaxy 
in the local universe) and the evidence for extinction in excess of the Galactic value, we consider a 
magnetar origin relatively unlikely. A long-duration GRB, however, can account for the majority of 
the observed properties of this source. GRB 090709A is distinguished from other long-duration GRBs 
primarily by the large amount of obscuration from its host galaxy (Apc,obs ^ 2 mag) . 

Subject headings: gamma rays: bursts - stars: neutron 



1. INTRODUCTION 

A variety of astrophysical sources are capable of pro- 
ducing short, intense flashes of high-energy emission de- 
tectable by the current generation of gamma-ray satel- 
lites. These sources span an incredible range of the ob- 
servable universe, from electri cal discharges associa ted 
with thunderstorms on Earth (jFishman et al.l fl994f ) to 
the deaths of the earliest known stars in the universe 
(jTanvir et al.ll2009t [Salvaterra et al.ll2009D . 

Gamma-ray bursts (GRBs) arc the most luminous 
class of these high-energy transients (£; so « 10 50 - 
10 52 ergs _1 ). At least two di stinct progenitor system s 
are thought to produce GRBs (|Kouveliotou et al.lll993h . 
It is now widely accepted that most GRBs with du- 
ration tgo > 2 s arise as a byproduct of the core col- 
lapse of massive star s (hereafter referred to as "long- 
duration" GRBs; e.g.. lWooslev fe Bloomll2006l and refer- 
ences therein). The origin of "short-duration" GRBs is 
still a hotly debated top ic. They likely arise from an older 
stellar population fe.g-.lGehrels et al.ll2005l: iHiorth . et alJ 
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2006), possibly due to the merger of a neutron star- 
neutron sta r (NS-NS) or black hole-neutron star ( BH- 
NS) binary (jEichler et al.lll989t iNaravan et al.lll992f K 

Soft gamma-ray repeaters (SGRs), on the other hand, 
are distinguished from GRBs by repeated outbursts 
with isotropic energy release of E^.iso 5, 10 46 erg (e.g . , 
IWoods fc Thompsorj[200l lOf^200% IMeregh ctti 2008) . 
The di scovery of per i odic oscillations from bright SGR 
flares (|Mazets et al.l I19791 ). along with the measure- 
ment of a spindown in their periods (|Kouveliotou et "all 
Il998h . allows for an estimation of the magnetic field 
strength. Unlike typical rotation-powered radio pulsars, 
SGRs (as well as their counterparts discovered in quies- 
cence, the anomalous X-ray pulsars, or AXPs) are likely 
powered by their intense m agnetic fields (B > 10 14 G; 
iDuncan fc ThompsonlfT992l) . Together, SGRs and AXPs 
are now thought to comprise a single class of young, 
highly magnetized neutron stars ( "magnetars" ) . 

At least three discoveries have in recent years 
challenged this simple classification picture. First, 
both GRB 060614 and GRB 060505 had high-energy 
durations in excess of a few seconds. Yet de- 
spite being quite nearby (z i=s 0.1), neither exhib- 
ited evidenc e for an associated supernova to quite 
deep limits dGehrels et al.1 120061: iGal-Yam et~aTl 120061: 
Delia Valle et alJ 120061 : iFvnbo et all 120061 : lOfek et all 
2007f h Second, many GRBs that would have been 
classified as having short duration (t 90 < 2 s) by the 
less sensitive BATSE satellite exhibit soft (and often- 
times faint) X-ray tails extending hun dreds of seconds 
in time when observed by Swift (e.g., iBarthelmy et alJ 
120051* iPerlev et all I2009H ) - perceived duration is, af- 
ter all, at once a redshift-, sensitivity-, and bandpass- 
dependent quantity. Finally, the X-ray and optical af- 
terglow light curves of GRB 070610 (t 90 = 5 s) dis- 
played dram atic flares on time scales as short as sev- 
eral seconds (|Kasliwal et al.ll2008t[Stefanescu et alJl2008t 
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ICastro-Tirado et al.ll2008T) . Though almost certainly of 
Galactic origin, the source of this variability is still poorly 
understood. 

Here we present observations of GRB 090709A, a 
high-energy transient whose classification as a tradi- 
tional long-duration GR B has been called into question. 
iMarkwardt et a l. ( 2009) reported the detection of quasi- 
periodic oscillations (period P = 8.06 s) in the high- 
energy light curve, typical of the observed properties of 
SGRs and AXPs. However, unlike all currently known 
magnetars, the localization of GRB 090709A is incon- 
sistent with both the Galactic plane and any nearby 
galaxy. Through a detailed analysis of the high-energy 
light curve, the broadband afterglow, and the immedi- 
ate environment, we attempt to uncover the origin of 
GRB 090709A. 

Throughout this work we adopt a standard ACDM cos- 
mology (Q A = 0.73; tt m = 0.27; H = 71 kms" 1 Mpc" 1 ). 
All quoted uncertainties are la (i.e., 68%) confidence 
intervals unless otherwise stated. Spectral and tempo- 
ral power-law in dices are provide d using the convention 
/„ oc t~ a v-P (jSari et all fl998D . UT dates are used 
throughout this work. 

In the fina l stag es of preparing this manuscript, 
iDe Luca et all (|2009f) posted an analysis of X-ray (XMM- 
Newton and Swift-XRT) and gamma-ray (Swift-BAT 
and Integral-SPl/ ACS) observations of GRB090709A. 
These authors reach largely similar conclusions regard- 
ing the origin of GRB 090709A, although they favor a 
somewhat larger distance for the event. We attempt to 
highlight both the differences and similarities between 
our work in what follows. 

2. HIGH-ENERGY PROPERTIES 

GRB 0907 09A was detected by th e Burst Alert Tele- 
scope (B AT: [B arthclmv et al.ll2005a|) on-board the Swift 
satellite dGehrels et al 1)2004) at 7:38:34 UTC on 9 July 
2009 {t ; Mo rris et all (20091) . In Figure [1] we plot the 
15-350 keV BAT light curve, binned with Is time reso- 
lution, obtai n ed fol lowing the prescription described by 
iButler et all (pOTll . The emission is dominated by a 
broad peak beginning at to and lasting approximately 
100 s. However, there is evidence for low- level variabil- 
ity well before the trigger (t « to — 70s), and the light 
curve appears to rise again at t « t a + 400 s. Shortly 
thereafter the s pacecraft slewed away due to an observ- 
ing constraint (Sakam oto et al.l 120091 ). Including data 
from to — 66 s to to + 509 s, we measure tgo = 345 ± 64 s 
(15-350 keV). Over this interval we find the spectrum is 
best fit by a power law with an exponential cutoff, with 
a = 1.06 ±0.14 and E p = 299±^keV { X 2 = 38.71 
for 55 degrees of freedom, d.o.f.). The corresponding 15- 

350 keV fluence is / = 4.59i°'26 x 10~ 5 erg cm" 2 , making 
GRB 090709A one of the brightest Swift events detected 
to date. 

Superposed on the overall ~ 100 s duration rise and 
decline of the dominant emission component, the high- 
energy light curve exhibits relatively strong fluctuations. 
In the inset of Figure [IJ we plot the residual emission 
after subtracting a smoothed version of the 15-350 keV 
li ght curve (boxcar binned by 10 s). 

IMarkwardt et all (|2009f ) first reported a search for a 
periodic signal from GRB 090709A, claiming a detection 
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Fig. 1.— Swift-BAT 15-350 keV light curve of GRB090709A, 
referenced to the time of the BAT trigger, 7:38:34 UTC on 9 July 
2009. The high-energy emission is dominated by a ~ 100 s long 
peak beginning at to, although fainter emission is visible both be- 
fore and after this episode. The inset shows the residual emission 
after subtracting a smoothed version of the light curve. 

of an excess in the power density spectrum (PDS) at 
P = 8.06 s. To estimat e the significance of this excess, 
IMarkwardt et all (|2009[ ) normalize the entire PDS by the 
noise in the observed frequency range 0.2-0.6 Hz. After 
correcting for the number of frequency bins examined 
(~ 2000) and the estimated number of bursts for which 
such an analysis could be performed (~ 100), the au- 
thors conclude the observed peak is highly significant 
(null probability ~ 10~ 6 ). The detection of this ap- 
parent periodicity was subsequen tly confirmed on several 
additional high-energy satellites fGolenetskii et alJl2009b 
iGotz et alJl2009b lOhno et al.ll2009f )~ 

In the left panel of Figure [2] we plot the un-normalized 
PDS of GRB090709A. The 15-350 keV light curve from 
to to to + 100 s was binned with a time resolution of 
10 ms, and de-trended by subtracting a smoothed ver- 
sion of the light curve ( 10 s boxcar). The peak noted by 
IMarkwardt et al] ([2009) is clearly visible at a period of 
8.1s. It is crucial to note, however, that the interpreta- 
tion of the significance of this peak depends sensitively on 
the assumed noise properties and the underlying shape of 
the PDS. We therefore examine this issue here in greater 
detail. 

The simplest strategy to infer the statistical signif- 
icance of this feature is to assume that the PDS is 
dominated by white noise (i.e., independent of fre- 
quency). The significance (signal-to-noise ratio, SNR) 
can then be estimated by normalizing the PDS by the 
obs erved scatter in a region devoid of features. Follow- 
ing |MarkwardFeE2] (|2OO90 . we normalize the PDS of 
GRB 090709A with respect to the observed scatter in the 
range P = 3-6.5 s (one of the noisier regions of the PDS). 
The result is shown as the dashed line in the right panel 
of Figure [2] After correcting for ~ 10 5 trials (~ 1000 
period bins using our 10 ms light curve and considering 
< P < 10s for - 100 long-duration GRBs), we find 
that the observed peak at P = 8.1 s is significant at the 
12rr level. 

The PDSs of GRB prompt emission, however, are not 
featureless. Examining the left panel of Figure [2j three 
distinct regions can be defined. For the shortest periods 
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Fig. 2.— Left: Un-normalized PDS of GRB 090709A. The PDS was constructed using the 15-350 keV BAT light curve from t to 
to + 100 s. A smoothed (10 s boxcar) version of the light curve was first subtracted to remove the overall rise and decline trend. Three 
clear regions are visible. For P < Is, the PDS is dominated by Poisson statistics (i.e., white noise). For P > 10s, the effects of windowing 
and de-trending dominate the error in the PDS. The power-law slope in the intermediate regime (a = 2.5^ intrinsic to the GRB prompt 
emission but simulating the effects of red noise, is clearly visible. The red line shows our model for the underlying shape of the PDS derived 
using a non-windowed light curve absent de-trending. Noise models assuming a flat PDS (green line, f° noise) and a power-law PDS (blue 
line, / 2 5 noise) are also shown. Right: The normalized PDS of GRB090709A. The dashed line assumes a flat underlying spectrum (white 
noise), while the solid line accounts for the intrinsic fluctuations (red noise). Significance intervals assuming 10 s trials are indicated by the 
horizontal dashed blue lines. The significance of the observed feature at P = 8.1s drops dramatically after accounting for the underlying 
power-law spectrum. 



(P < Is), the PDS is relatively flat and dominated by 
Poisson statistical fluctuations (white noise). For large 
periods (P > 10 s), the PDS is dominated by the total 
duration of the analyzed light curve (the window func- 
tion, ~ 100 s for GRB 090709A) and the de-trending al- 
gorithm (~ 10 s smoothing for GRB090709A). 

For intermediat e per iods (1 < P < 10 s), 
iBeloborodov et all (|1998f > have shown that the PDSs 
of the longest (i 90 > 100 s) BATSE bursts are well- 
fit by a power law with index a « 5/3 (where power 
oc P a ). Since the average power-law index was found to 
be in close agreement with the Kolmogorov law, it has 
been suggested that the prom pt emission may be relate d 
to turbulence in the outflow (|Beloborodov et ail [2000). 
Subsequent analysis has indicated, however, that the ex- 
act value of the power-law slope can vary significantly 
from event to event, and is determined predominantly 
by the shape and duration of individ ual pulse episodes 
within the GRB (|Suzuki et al J 120021) . PDS measure- 
ments of brigh t Swift events have been performed (e.g., 
GRB080319B. lBloom et~aDl2009f ) and show results con- 
sistent with those from BATSE. 

Regardless of its origin, it is clear that we must include 
the underlying behavior of the PDS when evaluating the 
significance of individual features. The underlying slope 
is not "noise" in the sense that it is not caused by the 
limitations of the measuring instrument; the observed 
power-law PDS is an intrinsic property of the GRB itself. 
Nonetheless, we can remove this contribution much as we 
would eliminate red (a — 2) or pink (0 < a < 2) n oise 
caused by our measuring device (e.g., Vaughan 2005). 

To estimate the power-law slope, we construct a new 
PDS for GRB 090709A, including all of the available 15- 
350 keV BAT data and without any de-trending. This 
helps to remove the noise at P > 10 s and provides 
a longer lever arm to calculate the PDS slope. For 



GRB090709A, we then find a « 2.5. 

In the right panel of Figure [2 we re-normalize the 
PDS of GRB 090709A using the underlying spectrum de- 
scribed above (solid line). While the same shape is visible 
in the PDS, the significance of the peak has dropped dra- 
matically. After correcting for the number of trials, we 
find that the feature at P = 8.1s is significant at only 
the ~ 2cr level. 

As stated previously, the inferred significance de- 
pends sensitively on the assumed noise properties. The 
value we derive for GRB090709A {a = 2.5) is signifi- 
cantly steeper than the a verage BATSE value of ~ 5/3 
(|Beloborodov et all 119981 ). To investigate further, we 
have performed a comparable analysis on all BAT GRBs 
with a fluence at least 70% of the value derived for 
GRB090709A (to provide sufficient SNR) and a du- 
ration igo > 70s (to provide sufficient sensitivity to 
P w 10 s). Only 6 additional events meet this crite- 
ria: GRBs 041223, 061007, 080319B, 080607, 090201, 
and 090618. The results of this analysis are shown in 
Table CD and Figure O 

Each event in this sample exhibits a peak in the PDS 
with period ~ 10 s and significance greater than 2a when 
normalizing by a flat PDS. However, the significance of 
these features drops dramatically when normalizing by a 
power-law spectrum with index a ks 2-2.5. 

The apparent periodicity of GRB 090709A was also 
reported by several additional gamma-ray satellites 
(iGolenetskii et alJ [200l iGotz et al.1 [20091 : (Ohno et all 
l2009f ). presumably performing a similar analysis to 
that of iMarkwardt et al.l (|2009D . We repeated our 
analysis on the publicly available Suzaku light curve 
of GRB 090709A, and find our results essentially un- 
changed: after normalizing by a power law with index 
a ks 2.5, the significance of the feature at P = 8.1s 
falls to ~ 2-3o\ This is not entirely surprising, though, 



4 



Cenko et al. 



TABLE 1 

Power Density Spectra of Bright Swift-BAT GRBs 



GRB 


Period a 


White SNR b 


PDS Slope 


Red SNR C 




(s) 








090709A 


8.1 


14.1 


-2.5 


5.0 


041223 


16.5 


2.0 


-2.6 


0.4 


061007 


7.6 


1.8 


-2.2 


1.2 


080319B 


8.3 


4.6 


-2.3 


0.8 


080607 


7.7 


7.3 


-2.2 


2.7 


090201 


8.3 


8.8 


-1.9 


3.5 


090618 


15.2 


3.1 


-2.2 


1.2 



a Period of the strongest feature observed in the PDS. 
b Single-trial signal-to-noise ratio of the strongest feature ob- 
served in the PDS assuming a flat (white) intrinsic spectrum 
across the entire frequency range (normalized to the region 
P = 3-6.5 s). 

c Single-trial signal-to-noise ratio of the strongest feature ob- 
served in the PDS assuming an intrinsic power-law spectrum 
(red noise). See text for details. 

as the PDS power-law index has been shown to be 
relatively independent o f bandpass for a given GRB 
(iBeloborodov et al.lll998h . 

iDe Luca et al.1 (( 2009) have performed a similar timing 
analysis of both the Swift-BAT and Integral- API/ ACS 
(jGotz et al.ll2009l ) light curves of GRB090709A. These 
authors confirm our primary result that the reported pe- 
riodicity in the BAT light curve at P = 8.1s is only 
significant at the < 3cr level. Interestingly, the signal at 
this frequency appears to be much weaker in the SPI- 
ACS data. 

To summarize, while the high-energy light curve 
of GRB 090709A shows tantalizing evidence of quasi- 
periodic oscillations with P = 8.1s, we do not take it 
to be a significant feature required by the data. It is cer- 
tainly the most compelling candidate for such behavior 
amongst the long-duration GRB population observed by 
Swift. However, when properly accounting for the shape 
of the underlying spectrum, we find the significance of 
the observed periodicity is not sufficiently large to be 
conclusive. Instead, we must look to the afterglow and 
environment of GRB 090709A to attempt to unveil its 
progenitor. 

3. FOLLOW-UP OBSERVATIONS 

The X-Ray Telescope (XRT; iBurrows et all 120051 ) on- 
board Swift began observing the field of GRB 090709A 
beginning at to + 74 s. A fading X-ray counter- 
part inside the BAT error circle was pr omptly iden- 
tified and reported to the community dMorris et al.l 
l2009t lOsborne et all [20091 : iRowlinson & Morris! 1200*91 1 
Cross-matching field sources detected by the XRT with 
cataloged near-infrar ed (NIR) positions from 2MASS 
(jSkrutskie et al.l l2~006) , we measure a localization for the 
afterglow of a = 19 h 19 m 42'/46, S = +60°43' 39'.'6, with a 
90% containment radius of 1'.'2 (J2000.0; see lButlerll2007l 
for details). 

We plot the X-ray afterglow of GRB 090709A in Fig- 
ure El which we have taken from the online compilation 
of N.R.B. 11 Like many GRBs in the Swift era, the XRT 
began observations of GRB 090709A while the prompt 

11 h ttp:/ /astro. berkeley.edu/~nat/swift; see Butler & Kocevski 
112007ft for details. 



gamma-ray emission was still ongoing, and the earliest 
X-ray observations extrapolate smoothly to the tail of 
the prompt emission. After 10 3 s, the X-ray light curve 
is well described by a single power-law decay with index 
a x = 1-38 ± 0.02 ( X 2 = 394.93 for 417 d.o.f.). However, 
there is a large gap in the X-ray coverage from t w 2 days 
to t « 10 days. A broken power law with 2 < i& < 10 days 
can also provide a reasonable fit. In this case, the ini- 
tial decay index is somewhat more shallow (ai t x ~ 1-2); 
the break time and post-break decay index are not well 
constrained. 

The X-ray spectrum at this epoch is adequately fit by 
a power law with index (3 X = 0.95±0.07 (x 2 = 178.58 for 
175 d.o.f.), although the inferred X-ray column (Ah = 
1.831q2i x 10 21 cm~ 2 at z = 0) is significantly in ex- 
cess of the Galactic value (Ah Gai = 6.6 x 10 20 cm~ 2 ; 
IKalberla et a!ll2005Tl . In the case of an intrinsic power- 
law spectrum, the requirement for extinction in excess 
of the Galactic value is significant at the 13.5cr level. 
Fits with only Galactic extinction require at least two 
blackbody components and still provide worse quality 
than an absorbed power law (kT\ — 0.32 ± 0.02 keV, 
kT 2 = 1.11 ± 0.05 keV; X 2 = 198.64 for 174 d.o.f.). These 
results are broadly consis tent with those reported by 
IRowlinson fc Morris! (|2009f ). 

Both the Peters Automate d Infrared Telescope 
(PAIRITEL; iBloom et al"1l2006bf ) an d the robotic Palo - 
mar 60 inch (1.5 m) telescope (P60; iCenko et all I2006D 
automatically responded to the Swift alert and began 
observations within 2 minutes of the trigger time. Ad- 
ditional observations were obtained at later times with 
the OMEGA2000 NIR camera on the 3.5 m telescope at 
Calar Alto Observatory (10 July 2009), with the OSIRIS 
instrument (Cepa et al. 2009, in preparation) attached 
to the 10.4 m Gran Telescopio Canarias (GTC) at the 
IAC's Observatorio del Roque de los Muchachos in La 
Palma (11 July 2009), and with NLRC 2 behind the laser 
guide star adaptive optics system (jWizinowich et all 
I2006D on the 10m Keck II telescope (17 July 2009). All 
images were processed using standard routines in the 
IRAF 12 environment. Photometric calibra tion was per- 
forme d relative to the USNO-B1.0 catalog (jMonet et al.l 
2003) in th e optic al, using the filter transformations of 
Jordi et all (120061) where app ropriate, and relative to 
2MASS (ISkrutskie et al.ll200l in the NIR. The results 
of this campaign are listed in Table and displayed in 
F igure El 

iMorgan et all <|2009D reported the PAIRITEL after- 
glow detection in the H and K s filters; subsequent 
re-analysis has revealed a detection in the J band as 
well (Figure 2} . Using early-time Ar s -band images and 
reference objects from the 2MASS catalog, we mea- 
sure a position for the afterglow of a = 19 h 19 m 42'.'64, 
5 = +60°43'39'.'3, with a 90% containment radius of 0'.'4 
(J2000.0). This location is at a distance of 1'.'4 from the 
center of the X-ray error circle, consistent within the 90% 
confidence regions. 

Combining the PAIRITEL NIR detections with the 
marginal r' and z' detections from P60, we find that 

12 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association for Research in 
Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Fig. 3.— The PDSs of other Swift long-duration GRBs. We have selected all BAT events with fluence > 70% that of GRB 090709A and 
duration tgo > 70s. As in the right panel of Figure [2] the dashed lines plot the normalized PDS assuming a fiat underlying spectrum using 
the noise properties in the range 3 < P < 6.5 s. All six events have a peak at P as 10 s with single-trial significance > 2a, though none is 
as strong as that derived for GRB090709A. After correcting for the underlying power-law spectrum (solid lines) and the number of trials, 
none of these features remain with significance > 3o\ 



the afterglow was extremely red. After correcting for 
Gala ctic extinction (E[B — V] = 0.09 mag; Schlcgcl et al. 
1998), we fit the nearly simultaneous r' z'JHK s photom- 
etry at t « 350 s to a power-law spectral energy distribu- 
tion (SED) model and find an extremely steep spectral 
index: j3 Q = 3.8 ± 0.3 ( X 2 = 4.16 for 3 d.o.f.). The 
temporal decay index in the optical is not very well con- 
strained by our observations; assuming the same value 
for all filters, we find ao ~ 0.9. 

We observed the field of GRB 090709A with the Very 
Large Array (VLA) 13 at a frequency of 8.46 GHz on 
July 11.42 and at 1.43 GHz on July 11.39. For both 
observations the array was in the "C" configuration. 
The flux-density scale was tied to 3C 286 or 3C 147 and 
the phase was measured by switching between the GRB 
and a nearby, bright, point-source calibrator. To maxi- 
mize sensitivity, the full VLA continuum bandwidth (100 
MHz) was recorded in two 50 MHz bands. Data reduc- 
tion was carried out following standard practice in the 
AIPS software package. No emission is detected at the 
NIR afterglow location at either frequency to 2a limits 
of A,(8.46GHz) < 70mJv, /„ (1.43 GHz) < 288 fiJy (see 
also lChandra fc Frai]|[2009] ). 

13 The Very Large Array is operated by the National Radio As- 
tronomy Observatory, a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, 
Inc. 



Finally, we have obtained a series of pre-GRB opti- 
cal images with the Palomar 48 inch (1.2 m) telescope 
from the Deep Sky catalog 14 . 32 images of the field of 
GRB 09 0709A were taken as pa rt of the Palomar-QUEST 
survey (|Diorgovski et al1l2008f ) over a time interval of a 
month in August 2005. The images were obtained us- 
ing an order-blocking filter with a cutoff blueward of 
A « 6100 A, which we calibrate relative to the Sloan Dig- 
ital Sky Survey i' band. A coadded image of these frames 
does not reveal any source at the location of the afterglow 
of GRB 090709A to a 3er limit of i' > 23.2 mag. 

4. A NEW SOFT GAMMA-RAY REPEATER? 

The detection of quasi-periodic oscillations with P = 
8.1s would strongly suggest an SGR origin. Unlike 
normal (rotation-powered) pulsars, the spin periods of 
known SGRs and AXPs fall in the range 2.0-11.8 s 15 . 
We therefore consider if the observed properties (outside 
the prompt high-energy emission) could themselves be 
consistent with an SGR flare. 

By definition SGRs undergo repeated high-energy 
outbursts. We have therefore searched the historical 
GRB catalogs of both the Interplaneta ry Network (IPN: 
lLaros et al.|[l997t iHurlev et alTll999aD and BATSE for 



14 See http://supernova.lbl.gov/~nugent/deepsky.html 

15 See' for example, the McGill SGK/AXP Oniine"Catal og at 
http: / /www. physics, mcgill.ca/ ~pulsar/magnetar/main.html 
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Fig. 4.— Left: PAIRITEL false-color (JHK a ) image of the field of GRB 090709A, taken at a mean epoch of t + 379 s. The afterglow 
is centered inside the white circle (5" radius), and is clearly much redder than the surrounding field stars. Stars A and B are marked for 
reference in all images. Center: Keck/NIRC2 K' image of the same field, approximately one week later. The black circle has a radius of 
5". Right: GTC/OSIRIS V image of the identical field on 11 July 2009. No sources are visible in the immediate environment (< 4") of the 
afterglow location to a 3d limiting magnitude of K' > 22.8 mag, i' > 25.5 mag. All images are oriented with north up and east to the left. 



an excess of localizations consistent with the position of 
GRB 090709A (both missions provided relatively coarse 
angular l ocalizations , requiring such a statistical anal- 
ysis; see lOfekl [20071 for details). Using the IPN cata- 
log from 12 November 1990 to 31 October 2005 16 , we 
found four IPN localizations consistent with the position 
of GRB 090709A. However, running the same search on 
random pointings with the same ecliptic latitude results 
in a median overlap of four GRBs. A similar analysis 
with short-duration (t 90 < 1 s) BATSE GRBs in the 
current catalog 17 also reveals no significant excess in the 
direction of GRB 090709A. 

SG R flares are broadly divided into three classes 
(e.g.. I Woods fc Thompsoi3l25M lMereghettill200l l. The 
rarest of these are the giant fl ares, with only a 
handful observed to date ( e.g., iMazets et al.l 119791: 

~d] or 



Hurley et alj Il999bt iFeroci et all 119991 : IMazets et ail 
1999tlHurlev et aTll2005l ). The high-energy light curves of 
giant flares are dominated by a short (~ Is), hard spike 
of gamma-rays. Thereafter, the emission softens and de- 
cays exponentially with quasi-periodic oscillations at the 
spin period of the underlying neutron star. Clearly the 
bright initial spike was not observed from GRB 090709A. 
However it is possible that this component could be non- 
isotropic and beamed away from our line of sight, as some 
bright SGR flares lacking the initia l spike have been ob- 
served (e.g.. lMereghetti et al.ll2009ft . 

Giant SGR flares have gamma-ray energy releases 
~ 10 44 -10 46 erg. Using our derived 15-350 keV fluence 
(§HJ, the observed system would fall at a distance ~ 0.1- 
1 Mpc (note that any beaming correction would only 
make the event more nearby, a possibility we discuss sub- 
sequently). Such distances imply an origin either in the 
halo of the Milky Way or in a nearby galaxy in the Local 
Group. 

Magnetars have relatively short characteristic ages 
(P/ P ~ 10 3 -10 5 yr), which would seem to rule out a halo 
origin if the neutron star formed directly from the core 
collapse of a massive star. If the progenitor system were 
associated with an older stellar population, for instance 



16 



See |http://www.ssl. berkeley.edu/ipn3/interpl a.html| 
version 17 December 2005. For reference, 

this IPN catalog version is also available at 



the accretion-induced colla pse (AIC) of a white dwarf 
([Canal fc S chatzmanl 119761 ) . this could perhaps explain 
the large offset from a site of recent star formation. How- 
ever, a relatively unphysical model would be required to 
account for the observed X-ray spectrum without excess 
extinction (see below). 

Giant magnetar flares from nearby galaxies are now 
believed to account for some fraction of the ob- 



Gaensler et al.l 120051 iCameron et al. 


2005; iTanvir et al. 


2005: Ofek et al.| 


20061 Ofek 2007 




Ofek et al.| 2008: 


Mazets et al.l |2008 


: iHurlev et all 12009 


; iChapman et al.l 



20091 ). In the case of GRB090709A, however, the lack 
of an obvious host candidate seems to disfavor an ex- 
tragalactic (but nearby; d < 1 Mpc) origin. Our pre- 
(Deep Sky) and post- (NIRC2) GRB limits at the lo- 
cation of GRB 090709A correspond to absolute magni- 
tudes of Mj/(AB) > -2 mag, M^-/(Vega) > -2 mag for 
d w 1 Mpc, sufficient to detect individual supergiant 
stars, let alone dwarf galaxies (modulo extinction). The 
post-GRB i' limits from the GTC provide even tighter 
constraints. 

Even with a large natal kick velocity (10 4 kms _1 ), the 
short lifetime limits the distance a magnetar can travel 
away from the host-galaxy disk (where it was presumably 
formed) to < 1 kpc. For d ~ 1 Mpc, this corresponds 
to an angular offset of ~ 3'. No Local Group galaxy 
is known within 10° of GRB090709A. The lack of any 
observed candidate host galaxy at this location seems 
to independently rule out a giant magnetar flare having 
extragalactic origin. 

Intermediate magnetar flares have durations of order 
1 s or longer and energy releases ~ 10 41 -10 43 erg. With 
the observed fluence of GRB 090709A, this corresponds 
to distances of ~ 4-40 kpc. The primary drawback 
with this scenario is the location of GRB 090709A with 
respect to the Galactic plane. With Galactic coordi- 
nates (l,b) = (91.8°, 20.1°), a total distance of > 4 kpc 
corresponds to a height of > 1 kpc above the Galactic 
plane. Of the 18 known and suspected magnetars, 
only two have Galactic latitudes 1 6] > 2°: SGR 0525-66 
(in the LMC: IMazets et all 119821: iKulkarni etai] 120031 : 
IKlose et al.l [2004J), and CXOU J010043.1-721134 (in 



|http://astro.caltech.edu/~eran/ GRB/IPN/NearbyGal/CatIPN .txt.verl71 22005 
11 See http:/ /gammaray.gsfc. nasa.gov/batse/grb/catalog/current 
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TABLE 2 

Optical/NIR Observations of GRB090709A 



Observation Date Time Since Burst a Telescope Filter Exposure Time Magnitude 15 Reference 





(2009 UT) 








fsl 






Jul 


9 07:40:22 


120.0 


PAIRITEL 


H 


23.4 


15.52 ± 0.16 


* 


Jul 


9 07:40:22 


120.0 


PAIRITEL 


K a 


23.4 


14.48 ± 0.21 


* 


Jul 


9 07:40:26 


142.5 


P60 


r' 


60.0 


21.80 ±0.51 


* 


Jul 


9 07:40:58 


177.0 


PAIRITEL 


H 


46.8 


> 15.90 


* 


Jul 


9 07:40:58 


177.0 


PAIRITEL 


K s 


46.8 


14.25 ± 0.14 


* 


Jul 


9 07:42:17 


307.5 


PAIRITEL 


H 


117.0 


16.20 ± 0.13 




Jul 


y U7:42:17 


307.5 


T~> A ID T 1 1 ' L/T 

FAIKI 1 -bL 


K s 


117.0 


15. UU ± U.15 


* 


In, 1 

JUl 


y u ( :4o: i i 


olo.z 


r^DU 


z 




1 n f^.A _i_ n Qn 
iy.o4 ± U.oU 


* 


hil 
J Ul 


Q n7-AD-^£ 


OOO. 


PA TTJTTTTT 
r\f\ir\l 1 1.1. 


1 




1 Q AQ _|_ fl 98 


* 


Jul 


Q n7-zL^-^^ 
y u < .^io. uo 


398.7 


P60 


r' 


180.0 


99 oa -\- n 49 


* 


Jul 


9 07:45:08 


484.0 


P60 


i' 


180.0 


> 21.62 


* 


Jul 


9 07:46:33 


569.4 


P60 


z' 


180.0 


> 19.92 


* 


Jul 


9 07:45:20 


680.0 


PAIRITEL 


H 


351.0 


16.88 ± 0.16 


* 


Jul 


9 07:45:20 


680.0 


PAIRITEL 


K s 


351.0 


15.89 ±0.21 


* 


Jul 


9 07:46:07 


1.637 x 10 3 


FTN 


i' 


520.0 


> 20.6 


1 


Jul 


9 07:41:16 


1.737 x 10 3 


FTN 


H 


720.0 


22.7 ±0.5 


1 


Jul 


9 08:06:34 


2.040 x 10 3 


P60 


r' 


720.0 


> 22.81 


* 


Jul 


9 08:11:06 


2.251 x 10 3 


P60 


i' 


600.0 


> 22.36 


* 


Jul 


9 08:13:11 


2.377 x 10 3 


P60 


z' 


600.0 


> 20.75 


* 


Jul 


9 08:20:17 


2.743 x 10 3 


P60 


</ 


480.0 


> 22.09 


* 


Jul 


9 09:30:00 


6.686 x 10 3 


Subaru / MOIRCS 


K 




18.2 


2 


Jul 


9 22:37:55 


5.479 x 10 4 


NOT / ALFOSC 


R 


1800.0 


> 24.0 


3 


Jul 


9 23:46:42 


6.488 x 10 4 


CAHA / OMEGA2000 


H 


3600.0 


> 20.73 


* 


Jul 


10 00:48:09 


6.777 x 10 4 


CAHA / OMEGA2000 


K s 


4140.0 


20.77 ± 0.39 


* 


Jul 


10 01:57:43 


6.801 x 10 4 


CAHA / OMEGA2000 


J 


3120.0 


> 21.74 


* 


Jul 


11 00:45:00 


1.480 x 10 5 


GTC / OSIRIS 


i' 


2100.0 


> 25.5 


* 


Jul 


17 10:06:16 


7.009 x 10 s 


Keck II / NIRC2 


K' 


1500.0 


> 22.81 


* 



References. — * = this work; 1 = IGuidorzi et all {2009); 2 = IAoki et"aTl (I2009D : 3 = lMalesani et al.l J2009D . 
a Time elapsed from the SW/i-BAT trigger, 07:38:34 UT on 9 J uly 2009 to the midpoint of the exposure. 

Magnitudes are provided in the AB system (Okc ik. Gunn 1983) for the g'r'i'z' filters and in the Vega system for the remaining 
filters. Errors quoted arc lcr photometric and instrumental errors summed in quadrature. For measurements with significance 
< 2o", 3<r upper limits arc provided. Galactic extinction (E(B — V) — 0.09 mag l lSchlegel eFal]|199a) has not been incorporated 
in the reported magnitudes. 



the SMC; lLamb et alJ l2002h 18 . For those SGRs in 
our Galaxy with known distances, the typical scale 
height is ~ 100 pc, as would be expected from their 
young ages. F urthermore, severa l are associated with 



SN remnants dKlose et all |2004t iGelfand fc Gaensleii 



200 

Vasisht 



Gaensler et all 
Gotthelf 



T9M iFahlman fc Gregory! I198jT 
or young stellar clus- 
iCorbel fc Eikenberrvl [200l 
even in the absence of a planar 

be 



ters (|Vrba et all 
iMuno et al.ll2006D . so 

birth, some remnant of recent star formation should 
evident nearby. 

lastly, we consider SGR flares with energy release 
< 10 41 erg. Such flares are quite common, though typi- 
cally observed with durations ~ 100 ms. The spectra of 
these flares are well described by an optically thin ther- 
mal bremsstrahlung model, with kT ~ 20-40 keV, unlike 
the observed properties of GRB 090709A. If we require a 
scale height of < 100 pc, this limits the total distance to 
GRB090709A to be < 300 pc. The energy release would 
therefore be < 5 x 10 38 erg, consistent with observations 
of low-level flaring activity from known SGRs and AXPs. 
However, the small distance would present additional dif- 
ficulties. Nearly all sufficiently well-localized magnetars 
exhibit quiescent X-ray emission at a luminosity of 10 33 - 
10 36 ergs~ 1 . At d < 300 pc, such emission should have 



18 The localization of an 
SGR0418±5729, is centered at 



has a radius of < 



additional SGR candidate, 
= 5.1°, b ut the error circle 
7° Ivan der Horst et alj|200Sft . 



been detecte d by X-ray surveys such as the ROSAT All- 
Sky Survey (jVoges et al.lll999l ). yet no cataloged X-ray 
sources are consistent with this position. 

Finally, we note that all Galactic models, even if as- 
sociated with an older population and not subject to 
the requirement of nearby star formation, would strug- 
gle to explain the excess extinction inferred from the X- 
ray spectrum and optical SED. While the intrinsic X-ray 
and optical spectra need not be power laws, we find no 
physically-motivated models capable of reproducing the 
observed X-ray spectrum and optical SED absent an ex- 
cess absorption component. 

5. MASSIVE STAR CORE COLLAPSE? 

Setting aside again the high-energy light curve, we con- 
sider the long-lived afterglow emission associated with 
GRB090709A. Radio afterglow emiss ion has be e n ob- 
served from two giant SGR flares dFrail et all 119991 : 
iCameron et all [20051: iGaensler et all I2005D . and fading 
X-ray counterparts are relatively comm on, even for in- 
termediate flares (e.g., Mercghetti 2008). However, the 
physical mechanism underlying this emission is thought 
to be quite different from the process powering long- 
duration GRB afterglows. 

In the case of magnetar flares, afterglow emission re- 
sults from a rapid injection of relativistic particles into 
the circumburst medium. The particles expand adi- 
abatically, and, consequently, the light curves decline 
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quite rapidly wi t h time (power- l aw indices a > 2.5; 
Frail et all 119991: ICheng fc Wand [2005 I Cameron et all 
20051: iGaensler et al.ll2005f) . For GRB afterglows, a rel- 
ativistic shock wave accelerates electrons in the circum- 
burst med ium, causin g them to emit synchrotron radia- 
tion (e.g.. IPiranll2005r ). As the shock expands outward, 
it continues to excite electrons further and further away 
from the explosion site, resulting in additional emission 
from larger rad ii and a slower fading rate (a 1-2; 
ISari et al.lll998h . 

A detailed comparison of well-sampled, broadband 
light curves may therefore be able to distinguish between 
these two scenarios. In the case of GRB 090709A, how- 
ever, the lack of optical/NIR and radio data make such 
a comparison difficult. We therefore attempt to answer 
a simpler question: is the observed emission consistent 
with standard GRB afterglow models? 

To do so, we r equire some knowle dge of the distance 
to GRB 090709A. lGrupe et all (|200l have demonstrated 
that the detection of X-ray absorption in excess of the 
Galactic value can be used to place a probabilistic upper 
limit on the distance to GRBs; the steep redshift depen- 
dence of soft X-ray absorption (~ (1 + z) 8 / 3 ) effectively 
precludes the observation of distant event s with large 
colum n densities. Using Equation 1 from iGrupe et al.1 
(l200l and our derived X-ray column excess, we find 
z < 3.5. We caution, however, that some well-known 
events are k nown to violate this redshift limit (e.g., 
GRB 080607; IProchaska et~alll2009h . 

A more strict li mit can be placed usin g the P60 and 
Faulkes Telescope (|Guidorzi et al.l 120091 ) r'-band after- 
glow detections. If we assume the Ly-a break in the 
GRB rest frame falls at the red edge of the observed 
i?-band filter, we conclude z < 4.5. 

An alternate approach is to use the observed BAT and 
XRT luminosity functions as a guide. The left panel of 
Figure [5] displays the X-ray afterglow flux (0.3-10 keV) 
at a common time of t = 1 day (observer frame) for all 
Swift GRBs to date with known redshifts. The X-ray 
afterglow of GRB 090709A is one of the brightest in the 
Swift era, well above any previous event with z > 4. 

Much like the X-ray flux, the BAT fluence from 
GRB 090709A is also one of the largest ever seen by 
Swift. In the right panel of Figure El we plot the rest 
frame l-10 4 keV isotropic energy release (E 7 ^ so ) for all 
BAT GRBs with known redshifts, along with K-corrected 
values for GRB090709A at z = 1, 2, 3, 4, and 5. 
To be consistent with the observed BAT distribution, 
GRB 090709A should fall at z < 3. More quantitatively , 
applying the GRB world model from lButler et all ([2009) , 
we find z < 2.5 (z < 1) at 99% (90%) confidence. 

An approximate lower limit can be derived using the 
lack of a host candidate at the location of the after- 
glow of GRB 090709A. W ith only a handful of exceptions 
(e.g. JCenko et al.ll2008l ). the observed host-offset distri- 
bution of long-duration GRBs is limite d to r < lOkpc 
(|Bloom et al.ll2002T : iFruchter et al.ll2006l ). Long-duration 
GRB hosts are typically blue, irregular gal axies (e.g., 
IFruchter et all 120061: iWainwright etal] 120071 ) with low 
metallicity (e.e.. iStanek et al.ll2006h and large s pecific 
star-formation rates (e.g.. lChristensen eTall200"l . If we 
take the low-metallicity dwarf Small Magellanic Cloud as 
a proxy for our GRB host (My ~ —16.6 mag), the non- 



detection in the i' band with the GTC requires z > 0.5. 

Combining these results, we believe the redshift of 
GRB090709A falls som ewhere in the r a nge .5 < z < 
3.0. For comparison, iDe Luca et al.l (2009) favor a 
somewhat larger distance (3.7 < z < 4.5; 90% con- 
fidence interval) based only on X-ray spectral fits as- 
suming an intrinsic power-law spectrum. The dis- 
tance limits derived here are not particularly constrain- 
ing. In what follows, we adopt a fiducial redshift of 
~ 1. The implied isotropic gamma-ray energy re- 
lease is then E ly i SO w 3.1 x 10 53 erg, and the observed 
host-galaxy limits correspond to Mj'(AB) > —18.6 mag 
and Mk (AB) > —19.5 mag. Both results are consis- 
tent with the observed long-duration G RB population 
(jButler et all 120071 : ISavaglio et al.ll2009D . Where appro- 
priate, we discuss the impact of varying the distance on 
the results. 

Clearly the X-ray afterglow provides the best sam- 
pling, so we begin our investigation there. After the 
first - 10 3 s, the X -ray decay is quite smooth. There 
is no evidence for dramatic variability on short time 
scales seen, for example, in the X-ray light curve 
of GRB 070610 dKasliwal et all 12001 iStefanescu etal] 
I2008t iCastro-Tirado et al.ll2008l ). Most strikingly, a de- 
tailed timing analysis of X-ray afterglow of both XRT 
and XMM-Newton observations of GRB090709A re- 
veals no evidence for perio dicity superposed on the 
overall power-law decline (jMirabal fc Gotthell 120091 : 
IDe Luca et al.l 120091 ). Such behavior is almost always 
seen in the decaying X-ray phase of SGR flares, and its 
absence is a challenge to any magnetar model. 

On the other hand, synchrotron afterglow theory pre- 
dicts a series of possible relationships between the spec- 
tral and temporal indices known as closure relations. For 
the X-ray afterglow of GRB 090709A, the observed val- 
ues agree well with theoretical predictions if the X-ray 
band falls below the cooling frequency (i.e., vx < v c ) 
and the shock wave is expanding into a constant-density 
circumburst medium (p oc r ): a = 3(3/2. It is relatively 
unusual to observe the cooling frequency above the X- 
ray bands at this e arly time, though n ot unprecedented 
(e.g., GRB 060418: iCenko et afll2009af ). We note, how- 
ever, that the case of the X-ray band falling above the 
cooling frequency (yx > v c ) is marginally acceptable, 
particularly if the shallower initial decay index from the 
broken power-law fit is used (light curve behavior above 
v c is independent of the circumburst medium density pro- 
file). The case of the X-ray bandpass falling below v c in a 
wind-like circumburst environment (p oc r~ 2 ) is strongly 
ruled out for both X-ray light curve fits. 

The observed optical spectral index (Po = 3.8) is much 
steeper than predicted by afterglow theory. Further- 
more, if we consider the last reported i^-band detection 
from Calar Alto at t « 0.8 day (when the X-ray light 
curve is dominated by afterglow light), the measured X- 
ray to optical spectral index is Pox ~ 0.26, well below 
the " dark" burst threshold (Pox < 0.5: Uakobsson et al.l 
120041 ). Both facts suggest a significant amount of optical 
extinction along the line of sight. 

In order to make the intrinsic K-band to X-ray flux 
ratio consistent with predicted values (Pox > 0.5), the 
observed -ftf-band flux must be suppressed by dust along 
the line of sight by at least A^^bs ^ 2 mag. However, if 
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Fig. 5. — Left: The 0.3-10 keV flux at a common time of t = 1 day (observer frame) for all Sure/t GRBs with known redshift. GRB 090709A 
is indicated with a dashed horizontal line. Right: The prompt gamma-ray (1-10 4 keV rest frame) energy release for all Swift GRBs with 
known redshift. GRB 090709A is shown for z = 1, 2, 3, 4, and 5. In both cases, GRB090709A is one of the brightest events in the Swift 
era, suggesting an approximate upper limit on the GRB redshift of z < 3. 



the synchrotron cooling frequency does indeed fall above 
the X-ray bandpass, the SED should be continuous from 
the X-rays through the optical; i.e., (3ox ~ 1- In this 
case, the required host extinction would rise sharply, to 
A K ,obs ~ 6 mag. 

We can estimate the rest frame V-band extinction 
(^4y,host) by assuming the optical SED remains largely 
unchanged from our measurement at t + 350 s (fto ~ 
3.8). If the optical to X-ray spectral index varies from 
0-5 < Pox ^$ 1-0, we find a host galaxy extinction of 
Av^ost « 4-8 mag (A^ ho st ~ 2-6 mag) for z « 1 (z w 3). 

Though most well-studied long-duration afterglow 
sightlines exhibit relatively s mall amount s of h o st ex- 
tinction ((Ay) pa 0.2 mag; iKann et alJ I2006L 120071 : 
ISchadv et al.l [2009), this result i s likely strongly b i- 
ased toward low-^4y sightlines (e.g.. lCenko et al.ll2009bl ). 
GRB afterglows with host Ay values as large as 
~ 3 mag have been reported in th e literature (e.g., 
iProchaskaet al.ll2009t iRol et al.ll2007h . and even larger 
values have be en inferred for GRBs without measure d 
redshifts (e.g., iTanvir et all 120081 : IPerlev et al1l2009afl . 
Like GRB 090709A, these highly obscured events of- 
ten masquerade as high-redshift candid ates at first 
(e.g., lAoki et al.l [20091 : iMorgan et all 12009ft . Large host 
columns are on occasion expected given the likely ori- 
gin of long-duration GRBs in giant molecular clouds 
in the disks of their host galaxies. Furthermore, the 
non-detection of the h ost galaxy at mm wavelengths 
(|Morgan fc Bow er 2009) is not particularly constraining, 
as most dark-burst host galaxies do not differ dramat- 
ically from the larger long-duration GR B population, 
suggesting a patchy dust distribution (e.g.. iBerger et all 
2003; IPerlev et al.ll2009aft . 

Examining a variety of different extinction laws over 
the redshift range 0.5 < z < 3, we find that a large selec- 
tive extinction coefficient (Ry ~ 4-15) is required to re- 
cover an intrinsic optical spectral slope of < (3o ^ 1 in 
nearly all cases. For comparison, the global val ue of Ry 
for the Milky Way (SMC) is 3.08 r2.93: IPeill992ft . A large 
value of Ry ( ~ 4) was inferred for the heavily obscured 
GRB 080607 (jProchaska et al.ll2009ft . similar to the lines 
of sight toward molecular clouds in the Milky Way, and 
may therefore be common for heav ily extinguished envi - 
ronments. The extinction law from[Calzetti e t al.l (|1994l ) 
can provide a reasonable fit across the entire redshift 



range of interest. To a large extent this is due to the 
featureless nature of this extinction law outside the far 
ultraviolet, which by default maintains the power-law 
shape of the observed optical SED. On the other hand, an 
SMC-like extinction law is only consistent with z > 1.5. 
Extinction laws with a strong 2175 A bump (Milky Way, 
LMC) also favor redshifts near our upper limit. 

We have also attempted to constrain the host galaxy 
column density (iVf^host) using the late-time (t > 10 3 s) 
XRT spectrum. Not surprisingly, the required column 
rises sharply as a function of redshift, from iVj^^ost = 
(3.5±0.5) x 10 21 cm" 2 at z « 0.5 to N Htlwst = (3.6±0.5) x 
10 22 cm~ 2 at z w 3. O ur X-ray spectral fits a re consistent 
with those derived in iDe Luca et al.l (|2009l ) given some- 
what lower preferred redshifts. At the smallest distances, 
the inferred dust-to-metals ratio ( Ay \host /Nh. host ~ 1 _ 2 
xl0~ 21 ) is in excess of the Galactic value. This result 
is inconsist ent with observations of previous GRB sight - 
lines (e.g.. IGalama fc Wijersl [20011 : 1 Watson et al.l i2006). 
where the dust-to-metals ratio is typically well below the 
Galactic value. At larger distances, ^l^host drops and 
Aff.host rises, lowering the inferred ratio to 6-17 xlO -23 
at z re 3. This result suggests GRB090709A may fall at 
the high end of our preferred distan ce scale. 

Assuming a Calzctti et al.1 (|1994ft extinction law with 
a large Ry can account for the optical SED, we have con- 
structed sample models for the X-ray, -R'-band, and radio 
afterglow em ission from GRB 090709A using the software 
described by I Yost et al.l (|2003ft . No single forward-shock 
model can account for all of the observed data; a sam- 
ple attempt is shown with dashed lines in Figure [H The 
parameters of the dashed-line model are relatively con- 
trived to provide an extremely small cooling frequency 
(v c < vo for t < 0.01 day), unlike what we inferred pre- 
viously based on the X-ray spectral and temporal indices. 
This is not entirely surprising, as we know the early X- 
ray emission is often dominated by othe r components 
besides the afterglow (jNousek et al.l[200rl ). In particu- 
lar, energy injection in this phase can significantly alter 
the subsequent evolution of the afterglow. 

If we only consider the emission at t > 0.1 day, when 
the afterglow emission from the forward shock appears to 
dominate the observed emission, we can easily find good 
models with v c > vx ■ One example is shown as the solid 
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Fig. 6. — The X-ray, X-band, and radio afterglow of GRB090709A. The dashed lines indicate a model that attempts to explain the 
entirety of the observations (the model has the following parameters: -Eke iso = 1-5 x 10 53 erg, 9 = 0.2, p = 2.01, e e = 0.5, eg = 0.2, z = 1, 
n = 40cm -3 , An,obs = 6 mag). While this model can account for most of the observed features, it implies an extremely small cooling 
frequency {y c < vo for t < 0.01 day). The solid lines shows a model with more reasonable physical parameters designed to explain the 
late-time (t > 0.1 day) data (the model has the following parameters: -Eke. iso = 1-5 X 10 53 erg, 9 = 0.1, p = 2.06, e e = 0.6, eg = 0.1, z = 1, 



n = 0.01 cm 



i-ff ,obs 



3 mag). The lack of optical and radio coverage at these times makes such models highly degenerate. 



line in Figure [6] Lower redshifts typically provide better 
fits to the bright X-ray flux. However, the lack of optical 
and radio data at these times leaves the basic physical 
parameters of such models largely unconstrained. 

6. CONCLUSION 

We now return again to our original question: what 
is the origin of GRB 090709A? A careful examination of 
the high-energy light curve indicates that the previously 
claimed detection of quasi-periodic oscillations was in our 
opinion overstated: with a proper accounting for the un- 
derlying spectrum, we detect the periodic signal with 
only marginal (~ 2tr) significance. Together with the lack 
of obvious environmental clues of recent star formation, 
it seems unlikely that GRB 090709A could be caused by 
an SGR-like outburst from a highly magnetized neutron 
star. 

Independent of the high-energy properties, we have 
shown that GRB 090709A was almost certainly a cosmo- 
logically distant event. Even if associated with an older 
stellar population and therefore free of the requirement 
to be nearby recent star formation, the X-ray and opti- 
cal extinction in excess of the Galactic value are difficult 
to reconcile with an origin in the Milky Way or nearby 
universe. The bright high-energy fluence and X-ray after- 
glow favor a more nearby event (z w 1), while the X-ray 
and optical extinction properties point to a more distant 
origin (z ss 3). Perhaps the true value falls somewhere 
in between. 

A long-duration GRB (and, by assumption, the core 
collapse of a massive star) , on the other hand, can natu- 



rally account for the majority of the observed properties 
of GRB 090709A: the large line-of-sight extinction, the 
lack of an obvious host galaxy, and the late-time after- 
glow decay. In many ways GRB 090709 A appears to be a 
relati vely close analog of GRB 080607 (|Prochaska et al.l 
2009): in particular the large host extinction (AK,obs ^ 
2 mag) and the large selective extinction (Ry > 4) are 
shared by both events. Such highly obscured GRBs may 
be relatively common, but are more likely to be discov- 
ered recently due to the prompt localization capabilities 
of Swift and rapid ground-based follow-up in the NIR. If 
so, it will be important to incorporate such optically dark 
events into systematic studies of t he GRB population to 
minimize observational bias (e.g.. iPerlev et al. 2009a). 

Subsequent observations could definitively resolve this 
issue: additional outbursts from this location would re- 
quire an SGR-like (i.e., non-destructive) progenitor. On 
the other hand, the detection of a faint host galaxy at the 
location of the afterglow (and in particular measurement 
of a substantial redshift) would provide even stronger ev- 
idence in support of our long-duration GRB hypothesis. 
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